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CHAPTER  I 


INTRODUCTION 

Electrically  small  antennas  have  been  the  subject  of  a  great 
deal  of  research  because  of  their  importance  to  communication  systems. 
For  many  reasons,  it  is  often  impractical  to  use  an  electrically  large 
antenna  which  has  characteristics  desirable  for  many  applications.  So 
the  necessary  electrically  small  antenna  often  becomes  the  limiting 
factor  in  the  capability  of  a  communication  system.  The  subject  of 
this  investigation  is  the  use  of  electrically  small  antennas  in  a 
direction  finding  system. 

The  problem  that  is  addressed  can  be  defined  by  enumerating  the 
physical  constraints  and  generally  desirable  characteristics  of  the 
system.  The  antenna  radiating  elements  or  array  are  to  be  electrically 
small;  that  is,  less  than  a  wavelength  in  extent  and  must  operate  on 
or  very  near  a  metallic  support  cylinder  of  modest  electric  size. 

These  geometrical  specifications  make  the  problem  a  conformal  array 
type.  The  system  is  to  perform  over  a  broad  frequency  band.  The 
antenna  elements  are  to  remain  mechanically  stationary  and  electric¬ 
ally  passive.  Signal  processing  is  to  be  simple;  that  is,  the  receiver 
is  capable  of  only  minimal  signal  manipulation.  Azimuthal  direction 
of  the  incoming  signal  is  to  be  found  without  knowing  the  frequency, 
phase,  or  polarization  of  the  signal.  The  constraints  do  not  completely 
specify  any  particular  direction  finding  system,  but  they  do  charac¬ 
terize  a  general  class  of  systems. 

The  investigation  described  herein  was  accomplished  with  the 
aid  of  a  digital  computer  and  two  approximate  techniques  of  solving 
electromagnetic  radiation  problems.  Chapter  II  briefly  presents  the 
concept  behind  the  piecewise-sinusoidal  reaction  matching  technique. 
Chapter  III  shows  the  development  of  a  hybrid  method  which  uses  modal 
solutions  for  an  infinite  cylinder.  The  research  was  conducted  by 
modeling  various  electrically  small  radiators  on  the  cylinder  and 
analyzing  the  structures  with  the  methods  of  Chapters  II  and  III  for 
the  purpose  of  investigating  the  direction  scheme  of  Chapter  IV. 

The  main  contribution  of  this  thesis  is  the  direction  finding 
scheme  discussed  in  Chapter  IV.  Thus,  Chapter  IV  develops  the 
proposed  direction  finding  system  by  discussing  the  consequences  of 
the  constraints  and  presenting  results  for  a  particular  radiator. 

Results  of  this  radiator  demonstrate  the  method  which  uses  electric¬ 
ally  small  elements  forming  an  array  around  the  cylinder  that  is 
capable  of  angle  of  arrival  determination  over  a  frequency  band  with 
minimal  signal  processing  when  the  frequency,  phase,  and  polariza¬ 
tion  of  the  signal  are  unknown. 
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CHAPTER  II 


SINUSOIDAL  REACTION  TECHNIQUE 


A.  Introduction 


Solutions  of  electromagnetic  radiation  problems  are,  except  for 
simple  geometries,  only  possible  through  approximations.  A  recent 
and  very  successful  approximation  method  involves  modeling  antennas 
or  scatterers  by  thin-wire  grids.  A  computer  program  for  thin-wire 
structures  in  a  homogeneous  conducting  medium  has  been  published  by 
Richmond  [1],  This  program  was  adopted  and  used  extensively  to 
analyze  the  various  antenna  structures  investigated.  The  program 
is  a  specialization  of  the  Method  of  Moments.  This  chapter  will 
briefly  discuss  the  concepts  behind  the  program  and  then 
present  some  results  for  simple  radiators  in  the  presence  of  a 
finite  length  cylinder. 

B.  Method  of  Moments 


Electromagnetic  radiation  problems  can  be  represented  by  integral 
equations  containing  field  terms  and  appropriate  inhomogeneous  source 
terms, either  actual  or  equivalent.  The  equations  seldom  lead  to 
tractable  closed-form  analytic  solutions.  The  Method  of  Moments, 
however,  is  a  general  approach  which  essentially  reduces  the  integral 
equation  to  a  system  of  linear  algebraic  equations  with  unknown 
coefficients  related  to  an  expansion  of  a  current  distribution. 

Consider  a  perfectly-conducting  body  with  surface  current  density 
j£Q  induced  on  it.  Applying  boundary  conditions  at  the  surface  of  the 
conductor  gives  for  the  total  electric  field 


+  -tan 


(2-1) 


where  E|an  is  the  tangential  component  of  the  incident  electric  field 
and  E|an  is  the  scattered  electric  field  radiated  by  the  equivalent 
current  density  J£q,  The  subscript  tan  will  be  dropped  with  the 
understanding  that  the  fields  are  the  tangential  components,  so  Eq. 
(2-1)  becomes 


(2-2) 
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To  relate  the  equivalent  current  density  Jgg  to  the  scattered 
electric  field  Ef,  a  linear  operator  L  is  defined 

UOjq)  -  -Es  (2-3) 


where  the  concepts  of  linear  spaces  and  operators  are  used.  Thus, 
Eqs.  (2-2)  and  (2-3)  give 


l<4q>  *  i1 


(2-4) 


For  a  typical  electromagnetic  radiation  problem  E.1  is  a  known 
source  and  J£q  is  the  unknown  response  to  be  determined.  L  will  be 
an  integral  operator  which  operates  on  a  unique  j£Q  to  give  the  known 
excitation.  That  is,  L-'  exists  such  that 

O^'f'd1)  .  (2-5) 

For  a  given  problem  the  domain  of  the  operator  must  be  defined 
by  determining  the  space  of  functions  J£q  on  which  it  operates  with 
the  resulting  functions  forming  the  range  of  the  operator.  The 
operator  L  can  be  considered  as  the  transformation  mapping  from  a 
subset  containing  J£q  to  one  containing  ET .  Besides  determining  the 
range  and  domain  of  the  operator,  it  is  necessary  to  formulate  an 
inner  product  <«J,E>,  which  is  a  scalar  called  reaction.  The  inner 
product  is  defined  such  that  it  forms  a  unitary  space.  The  following 


axioms  are  satisfied: 

<J,E>  ■  <E,J>  (2-6) 

<aj  +  00,  E>  =  a<0,E>  +  6<J,E>  (2-7) 

If  <J*,J>  >  0,  then  J  f  0  (2-8) 

If  '  <J*,J>  =  0,  then  J  =  0  (2-9) 


where  a  and  3  are  scalars  and  *  denotes  complex  conjugation.  The 
solution  to  Eq.  (2-5)  follows  from  the  four-step  procedure: 

1.  Expand  the  unknown  in  a  series  of  basis  functions, 
spanning  O^g  in  the  domain  of  L. 
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2.  Determine  a  suitable  inner  product  and  define  a  set  of 
weighting  functions. 

3.  Take  the  inner  product  and  thereby  form  the  matrix 
equation. 

4.  Solve  the  matrix  equation  for  the  unknowns. 


Once  the  unknown  current  coefficients  are  obtained,  the  radiation 
characteristics  such  as  far-field  patterns  and  impedance  are  easily 
obtained. 

A  specialization  of  this  general  Method  of  Moments  which  is 
particularly  convenient  for  electromagnetic  radiation  problems  is 
Galerkin's  Method.  Thus,  step  (1)  requires  the  expansion  of  the 
unknown  j£q  in  a  series  of  basis  functions  ,  0?,  ...  Jm  on  the 

surface  of  the  conducting  body  and  spanning  the  domain  of  L.  Jjtq 
can  be  represented  as 


l 

m 


Im  J 
m  — m 


(2-10) 


where  I  ' s  are  the  desired  unknown  complex  current  coefficients. 
Substituting  Eq.  (2-10)  into  Eq.  (2-4)  gives 


l(  [  imJ)  =  Ei  C 

“  m  ~m  — 
m 

which  is 

l  UJj  =  E1  (; 

t  m  Mn  — 
m 

by  the  linearity  of  L. 

*  Step  (2)  defines  a  set  of  weighting  functions  W ^ ,  W^,  •••  Hn 
in  the  domain  of  L  and  forms  the  inner  product 


I 

m 


m 


<Hn» 


(2-13) 


Galerkin's  specialization  is  that  the  weighting  functions  are 
chosen  equal  to  the  expansion  functions,  that  is,  W_  =  ol  .  Thus, 
Eq.  (2-13)  becomes  1 


I 

m 


I  <J  ,  L(JJ> 
m  — n  — m 


-  <-n»  I1’ 


(2-14) 
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where  the  inner  products  are  identified  as  reactions. 

Step  (3)  calculates  the  various  inner  products  which  may  be 
represented  in  matrix  form: 


<J1,  L(J])> 

<Jr  l(J2)> 

•"  ’iv  L<4>' 

*1 

<02,  L(01)> 

•  999 

•"  L<i)> 

I  •  •  • 

lZ 

.• 

9 

= 

<4.* 

L 

l(«L)> 

— m  -m  _ 

<Jl ,  E *> 
:J2,  £1> 

mL,  £l> 
— m  — 

(2-15) 


This  may  be  written  in  compact  notation 


[Z]  (I)  -  (V) 


(2-16) 


where  the  elements  are  considered  as  generalized  impedances,  cjrrents 
and  voltages,  respectively. 

Step  (4)  involves  standard  matrix  inversion  techniques  to  solve 
for  the  current  column.  The  simple  notation  used  to  represent  this 
general  moment  method  formulation  makes  the  solution  appear  deceptively 
simple.  Each  of  the  mz  impedance  elements  in  Eq.  (2-15)  suggests  at 
least  two  integrations,  one  when  applying  the  integral  operator  and  the 
second  dictated  by  the  inner  product.  The  next  section  discusses  a 
method  which  is  particularly  efficient  when  applied  to  thin  wire 
solutions  that  can  be  used,  for  example,  in  wire-grid  modeling. 

C.  The  Piecev.’ise-Sinusoidal  Basis: 

Re act  i  o  n-MatcKi  n~ 

In  this  section,  the  method  [2]  of  choosing  piecewise- 
slnusoldal  basis  functions  to  obtain  the  elements  of  th»  .mpedance 
matrix  [Z]  using  the  reaction  concept  will  be  outlined.  The  reaction 
concept  is  also  called  reaction-matching  and  can  be  demonstrated  by 
examining  a  dipole  antenna  with. a  test  source  placed  at  the  center. 

For  the  linear  antenna  and  M1 ,  the  original  sources  may  be  replaced 
wi  th 


and 


n  x  H 


(2-17) 
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(2-18) 


on  the  closed  surface  of  the  dipole.  By  the  equivalence  principle 
j£p  and  Mjtq  on  the  surface  generate  the  fields  E_  and  H_  outside  the 
cnpole  surface.  The  dipole  is  perfectly  conducting  so  M£q  vanishes 
on  the  surface.  Next,  a  filamentary  electric  line  test  source  is 
placed  inside  the  dipole  coincident  with  its  axis.  This  test  source 
is  reacted  with  the  equivalent  surface  source  of  E_  and  H.  Using  the 
term  "react"  implies  the  application  of  the  reciprocity  theorem. 


<a,b> 


[J(a)  •  E(b)  -  M(a)  *  H(b)]  d  V&  .  (2-19) 


Here  J(a)  and  M(a)  are  the  test  source  distributions  generating  E(a) 
and  Hla)  in  the  presence  of  b,  the  equivalent  sources.  Reac¬ 
tion  then  is  a  measure  of  the  coupling  between  the  two  sources  a  and 
b.  The  equivalent  surface  current  density  Jeq  can  be  represented  by 
an  equivalent  filamentary  source  parallel  to  the  test  source  and 
displaced  by  the  radius  of  the  dipole. 


Assuming  a  piecewise-sinusoidal  current  distribution  on  the 
filamentary  test  source 


It  iz) 


_  sin  k  (h  -  j z | ) 
sin  k  n 


(2-20) 


where  the  z-oriented  source  has  length  2h,  z  is  measured  from  the 
coordinate  origin  and  source  center,  and  k  =  2tt/a.  The  current  on 
the  equivalent  source  is  expressed  as 


I  (z)  = 


sin  k  (h  -  z | ) 
sin  k  h 


(2-21) 


where  I]  is  the  desired  unknown  current  coefficient  caused  by  generator 
voltage  V-j .  By  the  reciprocity  theorem,  the  solution  is  found  from 
I-jZ  =  V-j  where  Z  is  the  reaction  between  the  test  source  with  known 
field  E.  and  H.  and  the  equivalent  source  Jrn  of  the  true  fields  E_ 
and  H.  L  w 

z  ■  T7"  17  4*  ds  (2'22) 

This  concept  of  a  test  source  with  a  piecewise-sinusoidal  current 
distribution  with  known  field  being  reacted  with  the  equivalent 
true  source  distribution  on  the  surface  of  a  radiator  is  expanded 
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to  generate  an  impedance  matrix  from  which  the  unknown  current 
coefficients  are  calculated.  Consider,  for  example,  a  circular 
loop  antenna  which  may  be  modeled  by  a  polygon  of  thin  wires  forming 
a  tubular  surface.  A  filamentary  dipole  is  placed  in  the  center  of 
and  coincident  with  two  segments  of  the  polygon.  Reciprocity  is 
enforced  successively  between  the  test  dipole  and  each  of  the 
tubular-surface  dipoles  in  turn,  thus  generating  one  row  of  the 
impedance  matrix.  The  test  source  is  then  moved  around  the  loop 
to  create  the  other  rows.  Relating  this  back  to  the  general  Method 
of  Moments,  it  is  seen  that  the  integral  operator  and  inner  product 
have  been  defined  by  the  reaction  integral  Eq.  (2-19).  The  expansion 
and  weighting  functions  are  the  piecewise-sinusoidal  functions  making 
this  a  form  of  Galerkin's  method.  Thiele  [3]  has  discussed  the 
advantages  of  this  reaction-matching  technique  in  terms  of  accuracy 
and  rate  of  convergence  over  other  methods. 

D.  Calculation  of  Antenna  Characteristics 

Once  the  current  is  known  on  the  perfectly-conducting  body,  that 
is,  the  antenna  structure,  calculation  of  other  antenna  characteristics 
follow  directly.  The  current  distribution  for  the  piecewise-sinusoidal 
functions  will  simply  be  the  superposition  of  all  the  equivalent  tubular 
dipole  currents  defined  at  various  points  along  the  antenna.  Once  the 
terminal  currents  are  known,  the  terminal  impedance  follows  from 


zt  --  v!t 


(2-23) 


where  is  the  terminal  voltage.  Lumped  loading  may  be  included 
by  adding  the  load  terms  to  the  appropriate  elements  of  the 
impedance  matrix. 

To  calculate  the  far  field  anywhere  in  space  due  to  a  known 
current  is  straightforward.  For  a  filamentary  source 


E(r,e,<}>)  =  - 


I  (O  F(e,$)  dt  (2-24) 


along 
the  wire 
antenna 


where 


F(e,<j>)  =  EXP  jk[x‘  sin  8  cos  H  y'  sin  e  sin  $  +  z'  cose] 


and  (x' ,  y',  z')  denotes  the  source  coordinates  while  (©,❖)  denotes 
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the  observation  point  in  the  far  field.  The  above  equation  must  be 
applied  at  points  over  the  entire  antenna  structure  and  the  results 
superimposed  to  give  the  total  far-field  Ee(r,e,<f>)  and  (r ,<*» ) . 


The  power  gain  is  determined  by 


G(M) 


r  M2*  U2  3rg, 

30  lltl2  Rt 


(2-25) 


where  is  the  real  part  of  the  antenna  terminal  impedance. 

Determination  of  the  radiation  efficiency  requires  calculating 
the  input  and  radiated  powers.  The  time  average  power  input  is 


(2-26) 


The  radiated  power  is 


Pr 


llJZ  R 


(2-27) 


where  Rp  is  the  radiation  resistance  which,  for  the  lossless  perfectly- 
conducting  case,  is  equal  to  Rt.  The  radiation  efficiency  is  then 
defined  as 


(2-28) 


E.  Thin-Mire  Modeling 

The  concept  discussed  in  this  section  is  the  application  of  the 
sinusoidal  reaction  technique  to  not  only  the  antenna  but  also  to 
the  metallic  environment  or  support  structure  of  the  antenna.  The 
surfaces  of  metal  bodies  are  modeled  by  a  system  or  grid  of  wire 
segments.  The  resulting  generalized  polygon  obtained  with  this 
modeling  technique  essentially  has  the  capability  of  representing  a 
surface  of  arbitrary  shape.  The  only  practical  limitation  of  the 
wire-grid  modeling  procedure  is  the  size  of  the  impedance  matrix  [Z] 
which  must  be  solved.  The  matrix  size  is  directly  related  to  the 
electrical  size  of  the  surface  area  that  must  be  modeled. 
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Two  important  examples  using  this  technique  will  now  be  presented, 
the  results  of  which  will  be  referred  to  throughout  this  text.  Con¬ 
sider  a  perfectly-conducting  circular  cylinder  of  finite  length  with 
first  a  radial  monopole  attached  and  second,  a  dipole  element  located 
near  the  cylinder  and  oriented  coaxial  to  it.  The  radiators  are 
modeled  with  consecutive  straight  wire  segments  while  the  cylinder  is 
modeled  with  connected  hexagons.  Figure  2-1  shows  the  configuration 
with  dimensions  and  orientations.  This  investigation  is  concerned 
with  results  over  a  frequency  range.  To  accomplish  this,  the  physical 
dimensions  of  the  model  are  held  constant,  while  the  frequency  and, 
correspondingly,  the  electrical  dimensions  are  varied.  The  computer 
program  accepts  as  input  the  locations  of  all  points  which  are  the 
vertices  of  the  segments  formed  by  connecting  the  points,  the 
conductivity  and  diameter  of  the  wires,  and  the  frequency  of  opera¬ 
tion. 


The  program  first  sorts  the  segments  and  forms  the  needed  number 
of  modes.  Next,  the  impedance  matrix  is  calculated  by  assuming 
piecewise-sinusoidal  currents  as  test  sources  inside  the  wires  and 
as  the  basis  functions  on  the  wire  surfaces  and  then  enforces  the 
reaction  relationships  between  them.  The  symmetric  properties  of 
the  matrix  are  used  to  reduce  storage  and  running  time.  The  matrix 
is  inverted  and  the  simultaneous  algebraic  equations  are  solved  for 
the  unknown  complex  current  coefficients.  As  the  previous  section 
demonstrated,  the  desired  antenna  characteristics  follow  directly 
from  this  point.  Figures  2-2  and  2-3  show  the  far-field  patterns 
decomposed  into  phi  and  theta  components  for  a  representative 
frequency  in  the  two  principal  planes  for  first  the  radial  monopole 
and  then  the  axial  dipole. 
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Fig.  2-1..  Wire  grid  model  of  circular  cylinder  with  a 
radial  monopole  or  an  axial  dipole.  The 
orientation  and  coordinate  system  is  the  one 
used  for  the  patterns.  The  dimensions  are 
in  centimeters. 
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Fig.  2-2(a).  Phi  component  of  the  far  electric  field 
pattern  in  the  horizontal  plane  for  the 
radial  monopole  of  Fig.  2-1  located  on  the 
finite  length  cylinder. 
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Fig.  2-2(b).  Theta  component  of  the  far  electric  field 
pattern  in  the  horizontal  plane  for  the 
radial  monopole  of  Fig.  2-1  located  on  the 
finite  length  cylinder. 
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Fig.  2-2(c).  Theta  component  of  the  far  electric  field 
pattern  in  the  vertical  plane  for  the 
radial  monopole  of  Fig.  2-1  located  on  the 
finite  length  cylinder. 
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Fig.  2-3(a).  Theta  component  of  the  far  electric  field 
pattern  in  the  horizontal  plane  for  the 
axial  dipole  of  Fig.  2-1  located  on  the 
finite  length  cylinder. 


CHAPTER  III 


RADIATING  ELEMENTS  ON  CIRCULAR  CYLINDERS 
OF  INFINITE  LENGTH 


A.  Introduction 


The  previous  chapter  presented  a  method  of  solving  electro¬ 
magnetic  radiation  problems  by  wire-grid  modeling  of  the  radiating 
structure.  For  practical  reasons,  this  method  is  limited  to 
structures  which  are  not  electrically  large.  The  structure  is 
modeled  with  wire  segments  no  longer  than  A/4  and  surfaces  with 
grids  not  greater  than  about  .2x  by  .2A.  The  number  of  segments 
and  grids  determines  the  number  of  modes  and,  correspondingly,  the 
size  of  the  impedance  matrix  and  number  of  simultaneous  equations 
to  solve.  The  computer  storage  and  running  time  are  approximately 
proportional  to  the  third  power  of  the  number  of  equations. 
Electrically  large  structures  require  numerous  modes  for  accurate 
modeling.  Thus  computer  storage  limitations  and  running  time 
expense  will  limit  the  size  of  the  structure  which  may  be  analyzed 
by  this  method.  This  constraint  provides  the  motivation  for 
considering  radiating  elements  on  infinitely  long,  perfectly- 
conducting  circular  cylinders  for  which  a  modal  solution  exists. 

The  modal  solution,  when  combined  with  data  obtained  by  the  method 
of  Chapter  II  for  the  current  distribution  on  the  radiating 
element  itsel f,  provides  an  accurate  and  computationally  efficient 
means  for  calculating  far-field  patterns  of  radiating  elements  on 
an  electrically  long  cylinder. 

The  motivation  for  requiring  a  method  to  examine  the  patterns  of 
small  antennas  on  long  cylinders  will  become  apparent  in  the  next 
chapter.  However,  generally  speaking,  the  length  of  the  cylinder 
does  have  an  effect  on  the  elevation  plane  pattern  and  this  effect 
will  have  significance  in  Chapter  IV. 

B.  Modal  Solutions  for  Infinite 

Circular  Cylinders'  ' 

The  modal  solutions  for  the  total  field  at  any  point  in  space 
for  a  plane-wave  incident  on  a  perfectly-conducting  infinitely  long 
circular  cylinder  are  presented  in  reference  [4]  by  Burnside.  The 
approach  is  similar  to  that  of  Sinclair  [5];  however,  it  is  based 
on  the  cylindrically  separable  fields  presented  by  Richmond  [6],  The 
solutions  are  divided  into  two  orthogonal  polarizations,  theta  and 
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phi.  Figure  3-1  illustrates  the  orientation  of  the  coordinate  system, 
antenna  elements  and  polarizations  of  the  incident  planewaves.  The 
phi  components  of  the  incident  fields  in  the  x-z  plane,  that  is  for 
=  0, are 


•s  -  jk(z  cos  0.-  +  x  sin  6 4 ) 

E'  _  •  » 

$  "  p  e 

.  „  jk(z  cos  0.  +  x  sin  e.) 

H]  =  e  e  1  1 

o 

ZQ  is  the  free  space  impedance  and  k  =  2n/x. 

The  z-component  of  the  incident  H-field  can  be  written  as 
.  sin  e.  jkz  cos  e.  jkx  sin  0. 

V = '  ~TT  e  e 


(3-1) 

(3-2) 


(3-3) 


where  the  Bessel  function  representation 

jkx  sin  e.  r 
e  1  =  l  j" 


n=-« 


Jn  (kp  sin  e ^ )  e' 


(3-4) 


may  be  used  to  give 


.  sin  e.  jkz  cos  r 
H* - z — “  e  1  i  j"  Jn(kp  sin  e,)  e-W 

**  n  n=-co 


(3-5) 


Using  cylindrically  separable  fields  from  Table  3-1,  the  phi 
component  of  the  phi  polarized  incident  electric  field  is 

.  jkz  cos  0.  ?  . 

E'  =  -j  e  1  l  jn  J'  (kp  sin  e  )  e“jn*  (3-6) 

♦*  n=—  n  1 

.  The  form  of  the  scattered  field  can  be  written 

„  jkz  cos  e. 


E  =  j  e 

A  A  J 


1  I  a  H'(2) 


an  H'^  (kP  sin  e.)  e’in*.  (3-7) 


n=-« 


where  the  an's  are  unknowns  which  may  be  evaluated  with  the  appropriate 
boundary  condition.  On  the  surface  of  the  perfectly-conducting  cylinder 
the  tangential  components  of  the  total  E-field  must  vanish.  That  is. 
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Infinite  length  circular  cylinder  with  location 
and  orientations  of  the  infinitesimal  radial 
monopole  and  axial  dipole. 


TABLE  3-1 


TE  AND  TM  FIELDS  SEPARABLE  IN  THE 
CYLINDRICAL  COORDINATE  SYSTEM 

The  harmonic  electromagnetic  fields  listed  below  satisfy  Maxwell’s 
equations  in  a  homogeneous  source-free  region. 


TE  Fields 


TM  Fields 


f- 


-  (jo»P/p)  R  p*  1 

Ep  = 

R'  P  V 

juu  R1  P  Z 

0 

0/s* 

R  P'  V 
R  P  Z 

R'  P  V 

Hp  = 

(jue/p) 

RP'  Z 

(1/p)  R  P*  V 

B2  R  P  Z 

hJ  = 

”z  - 

-jwe 

0 

R'  P  Z 

The  time  dependence  eJU)t  is 

understood. 

R  is  a 

function 

r  a  IUIH.UUM  »  y  un  ij  j  aim  l.  \  3  a  1  untuun  ui  1.  uniy*  rr  llilcb 

indicate  differentiation  with  respect  to  p ,  4>  or  2.  The  functions 
satisfy  the  following  differential  equations: 

p  d.  ip/.'l  +  (g2  p2  -  m2)  R  =  0 


P"  =  -  P 


Z"  =  -h  Z 


where 


B2  +  h2  = 

2 

w  ue  ,  and  B  and  h  are  constants. 

Some  solutions 

of  these  differential 

equations  are 

listed  1 

R  (p)  =  Jm(sp) 

P(<J> )  =  cos  n#» 

Z(z)  = 

cos  hz 

"m  «“>> 

sin  m4> 

sin  hz 

H»  ’  <*>> 

ejhz 

h‘2)  «x>) 

e-jm$ 

e-jhz 
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(3-8) 


Etan  =  °  =  ^  (P  =  9)  +  ^  (P  =  3)  • 

Substituting  Eq.  (3-6)  and  (3-7)  into  Eq.  (3-8),  and  solving  for  a  , 
gives  n 


=  .n  Jn  (ka  s1n  9i> 
n  '  (ka  sin  e.. ) 


(3-9) 


The  total  field  at  any  observation  point  is  now  specified  by  the 
above  relations  in  conjunction  with  Table  3-1.  Of  particular  interest 
to  problems  solved  in  this  investigation  are  the  following: 

jkz  cos  e.  T 

E  =  S. -  J  2n  jn+1  J  (kp  sin  e.) 

*p  kp  sin  e.  n=l  n  1 


J'  (ka  sin  e.)  m 

-^-rrn - - -  HI  }  (kp  sin  e.)  sin  n($  -  <t.) 

H'(2)  (ka  sin  e.)  n  1  1 

n  1  ..  i  ■} 


E  =  0  . 

4>z 


(3-10) 

(3-11) 


For  theta  polarized  incident  plane  waves: 

jkz  cos  e.  » 

Eep  =  'j  cos  ei  e  lQ  en  j  Jn  <kp  sin  9i) 

0  (ka  sin  e.)  m  L 

-  - - - h'v  '  (kp  sin  e.)  cos  n(<j>  -  $,) 

(^  sin  9.)  n 


(3-12) 


jkz  COS  6.0= 

E  =  -  sin  e  e  L  en  j  Jn  (kp  sin  9i} 


J_(ka  sin  0.) 

— rxt - - -  H^;  (kp  sin  e.)  cos  n(<<>  -  ) 

(ka  sin  e.)  "  1  #  , 


(3-13) 


where 


,  .fl.  =  =  0 

n  |2,  otherwise 


(3-14) 
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These  components  of  the  electric  field  will  be  the  ones  needed 
to  calculate  the  far  fields  radiated  by  infinitesimal  radial  mono- 
poles  and  axial  dipoles  located  near  the  infinite  cylinder.  The 
antennas  are  very  short  in  terms  of  wavelengths  so  that  the  open 
circuit  voltages  they  receive  will  cause  negligible  distortion  to 
the  field  which  would  exist  in  their  absence.  The  open  circuit 
voltages  are  equal  to  the  product  of  the  effective  heights  of 
the  radiators  times  the  component  of  electric  field  intensity 
parallel  to  their  axes.  That  is. 


iJmonopole  =  EeP  <p  =  b>  +  ♦  E*p  (e  =  •>)]  (3-15) 

dipole  =  £[e  Eez  (p  =  b)  *  ♦  E*z  (p  =  b)l  (3-16) 


where  t  is  the  effective  length  of  the  monopole  or  dipole  without 
the  cylinder  present.  Here  the  antennas  are  located  at  <j>  =  0,  z  =  0 
and  the  phase  is  referred  to  the  center  of  the  cylinder. 


For  the  case  when  the  antenna  elements  are  transmitting,  the 
electric  field  is  given  by 


j  z0  1  e' 


£(r*Vi )  =  — Vx? —  MW 


(3-17) 


1  is  a  complex  current  excitation.  The  next  section  shows  how 
of  differeJt°shapes  ^  USed  t0  find  the  far  field  °f  larger  radiators 

6.  A  Hybrid  Method 

orevw  c^Ssions  lor  !he  distant  fields  radiated, given  in  the 
section, are  based  on  infinitesimal  elements  of  electric 
J0  determine  the  far  fields  of  an  arbitrary  antenna 
the  pre^ence  of  the  ^finite  cylinder,  integration  over 

done^ner  r^iaPK  ure  curnrnts  must  be  performed.  This  can  be 
c2^eS?s  weioIJpH  h  sSerpo51tlon  of  th*  fields  from  the  infinitesimal 
That  is  9hted  by  tbe  appropriate  complex  current  coefficient. 


j  Z, 


kr 


:  — r  xr  J,  ‘n  **  <91  zn> 

(3-18) 
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where  In  is  the  current  of  the  nth  element  which  has  the  appropriate 
effective  height  Dn  and  whose  position  is  located  by  $n,  Zn. 

In  order  to  find  the  In's  necessary  for  this  calculation,  the 
methods  of  Chapter  II  may  be  employed.  Chapter  II  describes  how  a 
discrete  current  distribution  on  an  antenna  element  located  on  a 
finite  cylinder  is  calculated  using  wire-grid  modeling.  The  assump¬ 
tion  that  the  currents  induced  on  the  antenna  elements  for  a  particu¬ 
lar  excitation  will  be  almost  identical  for  finite  and  infinite  cylinder 
cases  is  made.  This  is  especially  true  for  cases  where  the  radiating 
elements  are  resistively  loaded  as  those  which  are  investigated  in 
Chapter  IV.  So,  the  current  distribution  calculated  by  the  methods  of 
Chapter  II  can  be  used  to  find  the  complex  weighting  coefficients  I„. 
of  Eq.  (3-18).  '  n 

This  technique  is  called  a  hybrid  method  since  it  combines  the 
thin  wire  methods  of  the  previous  chapter  with  the  modal  solution 
of  this  chapter.  Thus,  the  arbitrary  radiator  is  divided  into 
segments  oriented  as  radial  monopoles  or  axial  dipoles  and  of  extent 
no  greater  than  .lx.  Currents  at  the  endpoints  are  found  using  the 
thin-wire  computer  program.  The  complex  current  coefficient  for  any 
given  segment  would  be  the  average  of  its  endpoint  currents  multiplied 
by  its  length  in  wavelengths  normalized  by  a  normalizing  factor.  The 
far  field  will  be  the  sum  of  all  the  fields  radiated  by  each  segment. 

D.  Application  of  the  Hybrid  Method 

The  hybrid  method  discussed  in  the  previous  section  provides  an 
efficient  means  for  calculating  the  far-field  patterns  of  radiators 
located  on  or  near  an  infinite  cylinder.  At  the  end  of  Chapter  II, 
the  results  of  the  sinusoidal  reaction  technique  applied  to  a  radial 
monopole  and  axial  dipole  on  a  finite  cyliner  were  presented.  The 
finite  length  of  the  cylinder  has  very  pronounced  effects  on  the 
far-field  patterns.  Computer  capability  makes  extending  the  cylinder 
length  by  wire-grid  modeling  impractical.  Figure  2-1  shows  the  cylinder 
length  that  was  modeled  to  be  35  cm  which  is  only  0.875X  at  750  fliz. 

To  find  the  patterns  of  radiators  on  longer  cylinders,  the  modal 
solution  of  this  chapter  could  be  applied.  The  infinite  cylinder 
will  be  a  good  approximation  to  finite  cylinders  of  a  few  wavelengths 
since  the  radiators  considered  are  electrically  quite  small. 

As  a  comparison  of  methods,  to  show  the  effects  of  the  finite 
cylinder,  and  to  provide  a  reference  for  subsequent  discussions,  the 
radiators  of  Fig.  2-1  are  analyzed  on  an  infinite  cylinder.  The 
equations  of  the  far-field  expressions  are  programmed  along  with 
the  means  for  computing  the  complex  weighting  factors  using  the 
currents  calculated  for  the  radiators  in  Chapter  II. 

Figures  3-2  and  3-3  show  the  far-field  patterns  decomposed  into 
phi  and  theta  polarizations  in  the  two  principal  planes  of  the  radial 
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monopole  and  the  axial  dipole,  respectively.  Comparing  these  results 
with  those  on  a  finite  cylinder.  Figs.  2-2  and  2-3  give  the  following 
insights.  Phi  components  of  electric  field  radiated  by  the  radial 
monopole  are  not  affected  by  the  shortness  of  the  cylinder.  Theta 
components  of  the  field  lose  their  symmetry  and  front- to-back  ratio 
in  the  elevation  plane.  The  axial  dipole  transmits  no  phi  component 
for  either  cylinder.  The  theta  components  are  not  greatly  different 
at  the  frequency  shown. 

In  the  following  chapter,  the  method  of  this  chapter  along  with 
the  method  of  Chapter  II  will  be  applied  to  investigate  a  new  direc¬ 
tion  finding  scheme. 
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Fig.  3-2(a) .  Phi  component  of  the  far  electric  field 
pattern  in  the  horizontal  plane  for  the 
radial  monopole  of  Fig.  2-1  located  on  the 
infinite  length  cylinder. 
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Fin.  3-2(b).  Theta  component  of  the  far  electric  field 
pattern  in  the  vertical  plane  for  the 
radial  monopole  of  Fig.  2-1  located  on  the 
Infinite  length  cylinder. 
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Fig.  3~3(b).  Theta  component  of  the  far  electric  field 
pattern  in  the  vertical  plane  for  the 
axial  dipole  of  Fiq.  2-1  located  on  the 
Infinite  length  cylinder. 


CHAPTER  IV 


PROPOSED  DIRECTION  FINDING  SYSTEM 


A.  Introduction 


One  objective  common  to  all  radar  and  direction  finding  (DF) 
systems  is  the  determination  of  where  the  target  is  relative  to  the 
system  location.  One  parameter  basic  to  this  determination  is  the 
azimuthal  angle  in  the  system's  reference  coordinate  system.  This 
chapter  discusses  a  DF  system  proposed  to  work  under  the  constraints 
enumerated  in  Chapter  I.  A  system  which  works  over  a  frequency 
range  is  demonstrated  and  its  limitations  discussed.  The  particu¬ 
lar  antenna  configuration  which  is  used  to  demonstrate  the  system 
is  one  of  several  which  were  investigated  and  is  not  an  optimal 
design  for  a  particular  DF  system.  It  was  chosen  simply  as  a 
vehicle  to  present  the  concept  of  the  proposed  DF  system. 

B.  Method  of  Angle -of  Arrival 
Determination 


Many  methods  of  angle -of -arrival  determination  are  used  in 
radar  and  DF  systems.  The  methods  depend  on  available  information 
about  the  incident  or  return  signal  and  the  capability  of  the  re¬ 
ceiving  system,  including  the  antenna,  to  process  this  information. 
Many  systems  use  highly  directional  antennas  which  electrically  or 
mechanically  scan  the  area  of  interest.  This  study  is  restricted  to 
electrically  small  antennas  which  are  stationary  and  necessarily  have 
broad  patterns  of  low  directivity.  Many  systems  use  frequency  and 
phase  information  analyzed  by  a  complex  processor  to  determine  angle 
of  arrival.  This  study  proposes  a  system  requiring  only  the  most 
elementary  process inq  which  does  not  need  to  know  the  frequency  or 
phase  of  the  incident  signal  within  a  broad  frequency  band.  Many 
systems  receive  only  certain  polarizations  of  incoming  signals, 
while  the  proposed  system  responds  to  any  arbitrarily  polarized 
signal . 

The  patterns  of  electrically  small  radiators  in  the  presence 
of  an  electrically  larger  cylinder  are  restricted  to  several  basic 
shapes.  This  result  is  because  any  practical  element  less  than  a 
quarter  wavelength  in  extent  has  the  same  pattern  shape  independent 
of  its  current  distribution.  In  the  azimuthal  plane,  two  basic 
shapes  result  for  the  two  orthogonal  polarizations;  theta  (vertical) 
polarization  and  phi  (horizontal)  polarization.  These  transmitting 
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patterns  will,  by  reciprocity,  result  in  received  or  induced  voltage 
bitterns  of  the  same  shape.  The  two  basic  received  voltage  patterns 
are  shown  in  Fig.  4-1.  Section  E  of  this  chapter  demonstrates  by  a 
specific  example  that  these  patterns  can  be  received  over  a  broad 
frequency  band.  Section  E  also  demonstrates  that  an  arbitrarily 
polarized  incoming  signal  can  be  decomposed  into  these  two  patterns. 

In  general,  determination  of  angle  of  arrival  using  received 
voltage  patterns,  such  as  in  Fig.  4-1,  could  be  accomplished  as 
follows:  Space  three  receiving  ports,  one  every  120°,  about  the 
cylinder.  Consider  first  the  theta  polarized  component  of  the 
received  voltage  pattern,  V-theta  ($).  This  pattern  has  one  simple 
maximum  and  is  symmetric  about  phi  equal  180°  line.  Three  such 
patterns  spaced  120°  apart  would  allow  unambiguous  angle-of-arrival 
determination  by  simple  port  voltage  comparison  processing.  These 
three  patterns  are  illustrated  in  Fig.  4-2(a).  Figure  4-2(b)  shows 
how  the  360°  of  azimuth  angle  can  be  divided  into  six  sections  by 
comparison  of  voltage  magnitude  levels.  This  figure  also  shows 
that  the  ratio  or  the  difference  of  largest  to  next  largest  port 
voltage  gives  a  unique  number  for  each  azimuth  angle  in  its  range. 
This  algorithm  would  be  quite  simple  to  implement.  It  is  not 
presented  as  the  best  or  optimal  scheme,  but  merely  to  demonstrate 
that  the  angle  of  arrival  may  be  uniquely  determined  from  purely 
theta  polarized  voltages  induced  on  the  three  ports. 
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The  same  method  can  be  applied  to  the  port  voltages  induced  by 
purely  phi  polarized  incident  signals,  V-phi  (<j>),  such  as  shown  in 
Fig.  4-1.  This  pattern  has  two  simple  maxima  and  is  symmetric 
about  phi  equal  90°  or  180°  line,  that  is,  the  pattern  repeats  every 
180°  instead  of  every  360°.  Dividing  the  azimuthal  plane  by  comparison 
of  voltage  magnitude  levels  gives  six  pairs  of  sections.  The  number 
obtained  by  forming  the  ratio  or  difference  will  be  the  same  for  two 
angles  180°  apart.  The  angle  of  arrival  obtained  from  the  purely  phi 
polarized  induced  port  voltages  is  either  correct  or  else  wrong  by  180°. 
This  may  not  be  a  serious  ambiguity  for  the  following  reasons:  In 
practical  systems  the  direction  finding  signal  will  rarely  be  purely 
theta  or  phi  polarized,  but  will  contain  components  of  both.  Which 
of  the  two  angles  obtained  from  the  phi  scheme  is  correct  can  be 
determined  by  the  angle  given  by  the  theta  scheme.  For  the  frequency 
ranges  considered  in  this  study,  theta  components  of  the  incoming 
signals  will  be  dominant  because  of  the  propagation  properties  of 
practical  media. 

Angle  of  arrival  can  be  unambiguously  determined  if  port  voltage 
patterns  illustrated  in  Fig.  4-1  can  be  obtained.  If  they  are  re¬ 
ceived  over  a  frequency  band,  this  method  determines  the  angle  of 
arrival  of  an  incident  signal  using  small  passive  antenna  elements 
without  knowing  the  frequency,  phase  or  polarization  of  the  in¬ 
coming  signal. 
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f*9.  4-2(a) .  Theta  component  of  the  received  voltage 
pattern  for  three  input  ports  spaced  120° 
apart  around  the  cylinder. 
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C.  Problems  With  Simple  Radiators 

The  previous  section  demonstrated  the  desirability  of  obtaining 
patterns  like  those  shown  in  fig.  4-1  over  as  broad  a  frequency  bar.d 
as  possible.  For  a  purely  theta  polarized  incident  signal,  any  antenna 
element  close  to  a  large  conducting  cylinder  that  has  a  vertical  arm  less 
than  a  quarter  wavelength  in  extent  (or  dipole-type  feed  less  than  half 
wavelength)  will  have  a  horizontal  plane  pattern  similar  to  V-theta 
(<J>)  of  Fig.  4-1.  Similarly,  for  a  purely  phi  polarized  incident 
signal,  any  antenna  element  close  to  a  large  cylinder  that  has  a 
horizontal  arm  less  than  a  quarter  wavelength  in  extent  will  have  a 
horizontal  plane  pattern  similar  to  V-phi  (<P)  of  Fig.  4-1.  This 
was  demonstrated  clearly  for  various  radiators  analyzed  by  the  method 
of  Chapter  II. 

These  facts  could  also  be  argued  as  follows:  Any  arbitrarily 
shaped  antenna  in  the  presence  of  the  cylinder  can  be  decomposed  by 
the  principle  of  superposition  into  radiators  with  only  vertical  or 
horizontal  components.  These  two  basic  elements  have  pattern  shapes 
that  were  illustrated  in  Chapters  II  and  III,  Figs.  2-2,  2-3,  3-2, 
and  3-3.  These  patterns  are  transmitting  far- field  patterns,  but  by 
reciprocity  the  receiving  port  voltage  plots  can  be  predicted.  The 
theta  component  of  the  vertical  dipole  field  pattern  in  the  horizon¬ 
tal  plane.  Fig.  2-3(a),  would  have  the  same  receiving  pattern  as 
V-theta  (0)  of  Fig.  4-1.  Similarly,  the  phi  component  of  the  radial 
monopole  field  pattern,  Fig.  2-2 (a )  would  have  the  same  receiving 
pattern  as  V-phi  ($)  of  Fig.  4-1. 

Superposition  of  these  basic  elements  and  their  corresponding 
field  patterns  could  be  used  to  model  any  arbitrary  vertical  and/or 
horizontal  radiating  element  with  the  resulting  horizontal  plane 
patterns  being  similar  to  V-theta  (<|>)  and  V-phi  ($)  of  Fig.  4-1. 

So  for  purely  phi  or  theta  polarized  incident  signals,  many  acceptable 
antenna  elements  exist  which  give  the  desired  patterns  of  Fig.  4-1. 

Consider  a  radiating  element  which  has  both  horizontal  and 
vertical  arms  (bent  monopole,  half-loop,  etc.).  If  the  incoming 
signal  is  purely  theta  or  phi  polarized,  then  by  applying  the 
appropriate  angle  of  arrival  determination  method,  this  type  of 
element  spaced  every  120°  would  suffice  for  the  DF  system.  If,  however, 
the .polarization  of  the  incident  signal  is  unknown  or  not  pure,  that 
is,  arbitrary,  then  the  scheme  breaks  down.  With  unknown  polarization, 
which  scheme  should  be  applied,  theta  or  phi,  cannot  be  determined. 

For  any  given  incoming  signal,  the  two  schemes  or  algorithms  give 
different  answers  for  the  angle  of  arrival.  If  the  polarization  is 
arbitrary,  somewhere  between  purely  theta  or  phi,  this  problem  is 
complicated  because  neither  scheme  would  give  the  correct  answer 
since  the  received  voltage  pattern  would  be  some  weighted  combination 
of  the  V-theta  ( <t> )  and  V-phi  (<p)  of  Fig.  4-1.  This  combined 
pattern  would  be  different  for  each  unique  polarization  between 
purely  phi  or  theta. 
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A  solution  to  this  problem  could  be  achieved  if  the  incoming 
signal  could  be  separated  into  its  phi  and  theta  components.  Ideally, 
two  elements  might  be  used  to  cover  all  possible  polarizations  if  one 
responded  only  to  theta  polarized  signals  and  the  other  to  only  phi 
polarization.  The  theta  sensitive  element  would  use  the  scheme 
described  by  Fig.  4-2.  The  phi  sensitive  element  would  use  a  similar 
but  shifted  algorithm  to  give  the  angle  of  arrival,  or  the  angle 
+180°,  as  described  in  Section  B. 


The  axial  or  vertical  dipole  of  Chapter  II,  Fig.  2-1,  offers  the 
desired  theta  sensitive  element  with  negligible  phi  sensitivity.  A 
search  for  a  radiator  in  the  presence  of  the  cylinder  which  is 
sensitive  to  only  phi  polarized  incident  signals  from  all  an¬ 
gles  poses  a  difficult  problem.  Any  element  containing  horizontal 
arms  near  the  cylinder  produces  theta  components  of  the  far  zone 
electric  field  which  for  some  elevation  angles  are  as  strong  as  the 
phi  components.  This  was  investigated  for  various  radiators,  including 
radial  monopoles,  circumferential  dipoles,  circumferential  slots,  and 
loops.  So,  in  the  presence  of  the  cylinder,  two  elements  which  are 
sensitive  to  the  two  major  polarizations  independently  are  unavailable. 

The  next  section  demonstrates  that  the  proper  combination  of 
two  elements  facilitates  the  cancellation  of  components  of  one 
polarization. 

D.  A  Method  for  Separating  Polarizations 

To  achieve  the  desired  patterns  of  Fig.  4-1  over  a  frequency 
band, and  thus  form  a  DF  system,  the  incoming  signal  must  be  decom¬ 
posed  into  the  two  orthogonal  polarizations.  For  the  reasons  indi¬ 
cated  in  the  previous  section,  this  is  not  possible  using  radiating 
elements  which  are  independently  sensitive  to  the  two  polarizations. 
However,  two  symmetric  elements  correctly  positioned  and  properly 
combined  do  provide  the  desired  mechanism  for  separating  the 
polarizations. 

Consider  a  pair  of  vertically  oriented  radiators  such  as  half¬ 
loops  mounted  on  the  cylinder  as  illustrated  in  Fig.  4-3.  One  of 
the  pair  is  mounted  with  the  input  port  at  the  bottom,  the  current 
being  driven  to  the  top.  The  other  is  mounted  as  closely  as  possible 
to  the  first  element  but  inverted  so  that  it  is  fed  at  the  top  and 
terminated  at  the  bottom.  Careful  study  of  the  field  patterns  or 
induced  port  voltage  patterns  shows  that  the  two  elements  have  the 
same  phi  component  patterns  which  are  in  phase.  Theta  component 
patterns  are  also  the  same,  at  least  in  the  azimuthal  plane,  but 
are  180°  out  of  phase.  Examination  of  the  element  orientations  and 
physical  insight  justifies  these  results.  The  pair  of  elements  has 
vertical  components,  producing  theta  fields,  which  are  driven  in 
opposite  directions,  180°  out  of  phase.  The  pair  of  elements  has 
radial  components,  producing  phi  fields,  which  are  driven  in  the 
same  direction  and  are  therefore  in  phase.  The  radial  components  are 
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Fig.  4-3.  Two  resistively  loaded  half-loop  radiating 
elements  symmetrically  located  on  the  finite 
cylinder.  The  dimensions  are  in  centimeters 
and  the  coordinate  systems  are  the  ones  used 
in  the  patterns  which  follow. 


displaced  from  each  other  in  the  2  direction,  but  this  does  not 
greatly  affect  their  phases,  since  phi  components  are  independent 
of  z. 

If  the  ports  of  these  pairs  of  elements  are  considered  as 
input  terminals,  then  incoming  signals  will  induce  port  voltages 
which  can  be  combined  to  separate  the  two  polarizations.  Adding 
the  two  port  voltages  will  result  in  cancellation  of  the  theta 
components  of  the  incoming  signal  since  they  are  equal  but  180° 
out  of  phase  on  the  two  ports.  The  in-phase  and  equal  phi  components 
will  add,  giving  twice  what  appears  on  each  port.  Similarly,  sub¬ 
tracting  the  two  port  voltages  will  cancel  the  phi  contributions,  and 
double  the  theta  component  of  one  port  taken  alone.  For  this  method 
to  be  effective,  the  elements  must  be  symmetric  so  that  their  pat¬ 
tern  magnitudes  are  equal.  They  must  be  placed  as  close  together 
as  is  practical  so  that  their  phase  centers  are  as  nearly  superimposed 
as  possible  over  the  frequency  range.  It  is  important  that  the  ports 
maintain  their  relative  phase  difference  over  the  frequency  range  and 
for  all  aspect  angles. 

The  next  section  will  demonstrate  how  this  method  works  for  some 
pairs  of  radiating  elements.  By  decomposing  the  incoming  signal  into 
its  two  basic  polarizations,  it  becomes  possible  to  determine  the 
angle  of  arrival.  The  process  would  involve  two  steps.  First, 
adding  the  two  port  voltages  of  the  three  pairs  of  elements  spaced 
around  the  cylinder  would  give  only  the  phi  components  of  the  incoming 
signal.  These  three  voltage  levels  would  be  compared  as  discussed 
in  Section  B  to  give  the  angle  of  arrival.  Second,  subtracting  the 
two  voltages  of  three  pairs  of  elements  loaves  only  theta  components 
of  the  incoming  signal.  These  would  be  compared  to  give  the  angle 
of  arrival.  By  this  method,  the  desired  voltaoe  patterns  of  Fig.  4-1 
can  be  achieved,  and  thus  the  angle  of  arrival  can  be  determined  for 
a  signal  of  unknown  frequency  and  arbitrary  polarization. 

E.  Results  of  the  Method 


Figure  4-3  shows  a  pair  of  half-loops  attached  to  a  cylinder. 
This  type  of  antenna  element  was  one  of  several  studied  in  the 
presence  of  an  electrically  larger  conducting  cylinder  and  was 
chosen  to  demonstrate  the  method  of  separation  of  polarization 
discussed  in  the  previous  section.  The  half-loops  are  resistively 
loaded  to  achieve  desired  input  impedance  characteristics  as  dis¬ 
cussed  in  the  next  section.  The  DF  system  is  proposed  to  work  over 
as  large  a  frequency  band  as  possible.  To  study  this,  the  physical 
dimensions  of  the  cylinder  and  half-loops  are  held  constant  and  its 
antenna  characteristics  are  analyzed  over  a  frequency  range. 

For  the  fixed  dimensions  given  in  Fig.  4-3,  the  half-loops 
exhibit  desired  characteristics  over  the  frequency  range  250  MHz 
to  750  MHz,  a  three-to-one  band.  The  exact  bandwidth  would  depend 


_  c yctem  specifications  such  as  minimum  antenna  efficiency, 
l  Allowable  mismatch,  signal-to-noise  specifications  of  the 
T-nVsiqnal,  desired  accuracy  of  the  angle-of-arri val  determina- 
’  \  . -roC.*ssing  capability  of  the  receiver,  etc.  This  study  does 

I  .  Vie  exact  specifications  of  these  kinds  and  only  considers 
--••'oral  characteristics  desirable  to  all  DF  systems.  The  half-loops 
/'^-ified  in  Fig.  4-3  are  not  an  optimal  or  final  design  element 
"'/a  particular  system.  They  do,  however,  demonstrate  the  success 
cf  the  concept  of  the  DF  system  discussed  in  this  chapter. 


Figure  4-4  shows  the  received  voltages  combined  as  a  function 
c'  azimuthal  angle  in  the  horizontal  plane  as  described  in  Section  D. 
?-ree  frequencies  are  included  to  demonstrate  how  the  patterns  behave 
cvr  a  band.  The  incoming  signal  is  45°  polarized,  that  is,  it  contains 
....  phi  and  theta  components.  This  is  accomplished  as  follows.  First, 
. t,i  induced  port  voltages  at  Port  1  and  Port  2  are  calculated  by  the 
-••v-ods  of  Chapter  II  for  a  theta  polarized  incident  wave  of  magnitude 
1.0  V/m.  The  induced  port  voltages  for  a  phi  polarized  incident  wave 
ire  then  calculated.  The  complex  voltages  obtained  for  the  two 
polarizations  are  then  added  at  Port  1  and  Port  2.  The  voltages  at 
Pert  1,  V](*),  and  at  Port  2,  V2 (<(>),  are  those  that  would  be  induced 
:  /  a  45°  polarized  incident  wave  of  magnitude  /2  V/m.  The  polariza¬ 
tions  nay  now  be  separated  by  the  method  of  Section  D.  Addition  of 
voltages  at  Port  1  to  those  at  Port  2  gives  V-phi  (<f>),  and  their 
difference  gives  V-theta  ($). 


The  above  processes  can  be  described  by  the  following  equations: 


^phi  Port  1  ^  +  ^theta  Port  1  ^  ~  V1  ^ 

Vphi  Port  2  ^  +  V theta  Port  2  ^  “  V2  ^ 

Vphi  (4)  =  V]  (*)  +  V2  (♦) 

Vtheta  (♦)  "  V,  (♦)  -  V2  (♦) 

Figure  4-4  shows  the  results  of  this  process  for  three  frequencies 
in  the  horizontal  plane.  The  patterns  have  the  same  shape  as  those  of 
the  desired  patterns  which  Section  B  discussed  as  being  capable  of 
angle-of-arrival  determination.  Note  that  the  theta  polarized  voltage 
patterns  all  have  the  desired  single  maximums  which  occur  at  the  same 
'/i-’uth  angle  for  the  frequency  range  studied.  These  patterns  are 
‘  ‘  "'ll;/  suited  to  the  determination  of  incidence  angle  by  the  method 
described.  The  phi  polarized  patterns  are  not  quite  as  constant  over 
the  frequency  range.  The  maxima  occur  over  a  range  of  15  degrees  of 
azimuth  angle.  This  reduces  the  accuracy  of  the  algorithm  used  to 
determine  the  angle  of  arrival  for  this  polarization  by  approximately 
b  degrees  or  less.  Depending  on  the  specification  of  a  particular 
this  may  present  a  problem  or  be  acceptable.  So,  in  the 
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Fig.  4-4(a).  Theta  component  of  received  voltage  pattern 

obtained  by  subtracting  voltage  at  Port  1  from 
Port  2.  Three  frequencies  are  shown.  The 
incident  wave  is  in  the  horizontal  plane.  The 
half-loop  elements  are  on  the  finite  cylinder. 
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horizontal  plane,  the  method  of  separation  of  polarization  applied 
to  the  pairs  of  half-loop  elements  permits  angie-of-arri val  determina¬ 
tion  over  a  three-to-one  frequency  range  for  incoming  signals  of  any 

polarization. 

Extension  of  the  concept  of  separation  of  polarization  for 
determination  of  azimuthal  angle  of  arrival  for  incoming  signals 
from  some  elevation  angle  other  than  the  horizon  presents  some 
interesting  difficulties.  Figures  4-5  and  4-6  show  what  happens 
when  the  same  method  is  applied  to  incoming  signals  from  elevation 
angles  higher  than  the  horizon.  Figure  4-5  shows  the  combined 
voltage  patterns  for  an  incoming  signal  10°  above  the  horizon.  The 
patterns  have  shifted  slightly  where  maxima  and  minima  occur. 

This  would  result  in  a  corresponding  decrease  in  the  accuracy  of 
determined  angle  of  arrival.  Figure  4-6  shows  that  for  an  elevation 
angle  of  20°,  the  patterns  have  changed  and  shifted  enough  that 
applying  the  algorithms  would  result  in  fairly  large  errors. 

The  reason  for  the  problems  as  elevation  angles  are  increased 
is  due  to  asymmetries  in  the  elevation  plane  patterns  of  the  half¬ 
loop  radiator.  This  is  most  clearly  demonstrated  by  looking  at  the 
transmitting  elevation  plane  field  patterns.  Specifically,  Fig. 

4-7  shows  the  theta  component  of  the  far  electric  field  in  the 
elevation  plane  for  a  single  half-loop  excited  from  bottom  to  top. 

These  patterns  were  calculated  using  the  method  of  Chapter  II,  the 
finite  cylinder  case.  As  illustrated,  over  the  frequency  range, 
the  patterns  are  not  symmetric  about  a  line  for  theta  equal  to  90° 
representing  the  horizon.  A  half- loop  excited  from  top  to  bottom 
would  have  the  patterns  of  Fig.  4-7  but  inverted  by  180°.  The 
asymmetry  about  the  horizon  causes  difficulties  for  elevation  angles 
when  applying  the  method  of  separation  of  polarization  discussed  in 
Section  D.  These  figures  show  transmitting  patterns  but  by  recipro¬ 
city  the  induced  port  voltages  are  proportional,  that  is,  they  have 
the  same  shape. 

Consider,  for  example,  the  value  of  E-theta  at  750  MHz  for  an 
elevation  angle  30°  produced  by  the  half-loop  excited  from  bottom  to 
top.  Figure  4-7(c)  shows  the  magnitude  to  be  approximately  -1  dB. 

An  element  excited  from  top  to  bottom  would  have  magnitude  approxi¬ 
mately  -9  dB  as  given  at  theta  equal  120°.  The  method  described  for 
separation  of  polarization  depends  on  these  values  being  equal  in 
magnitude  and  180°  out  of  phase.  The  magnitudes  for  the  30°  elevation 
angle  are  quite  different  and  cancellation  of  the  theta  components  of 
combined  voltages  would  be  far  from  complete.  It  is  apparent  then 
that  the  application  of  this  method  of  polarization  separation  for 
elevation  angles  much  higher  than  the  horizon  requires  elevation 
plane  patterns  which  are  symmetric  about  the  horizon. 

Figure  4-7  displays  another  characteristic  of  the  theta  compon¬ 
ents  oT  the  far  electric  field.  At  750  and  500  MHz,  some  deep  nulls 
appear  for  various  small  elevation  angles.  These  nulls  cast  some  doubt  on 
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Fig.  4-5(a).  Theta  component  of  received  voltage  pattern 

obtained  by  subtracting  voltage  at  Port  1  from 
Port  2.  Three  frequencies  are  shown.  The 
incident  wave  is  10°  above  horizontal  plane.  The 
half- loop  elements  are  on  the  finite  cylinder. 
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Fig.  4-5(b).  Phi  component  of  received  voltage  pattern 
obtained  by  adding  voltage  at  Port  1  and 
Port  2.  Three  frequencies  are  shown.  The 
incident  wave  is  10°  above  horizontal  plane 
The  half-loop  elements  are  on  the  finite 
cylinder. 
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Fig.  4-6(a).  Theta  component  of  received  voltage  pattern 

obtained  by  subtracting  voltage  at  Port  1  from 
Port  2.  Three  frequencies  are  shown.  The 
incident  wave  is  20°  above  horizontal  plane. 
The  half-loop  elements  are  on  the  finite 
cylinder. 
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Fig.  4-6 ( b ) .  Phi  component  of  received  voltage  pattern 
obtained  by  adding  voltage  at  Port  1  and 
Port  2.  Three  frequencies  are  shown.  The 
incident  wave  is  20°  above  horizontal  plane. 
The  half-loop  elements  are  on  the  finite 
cyl inder. 
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the  practical  value  of  the  method,  since  deep  pattern  nulls  can  be 
very  sensitive  to  environmental  effects  in  a  practical  system. 

Small  variation  of  the  null  angles  could  reduce  the  stability  of 
the  front-to-back  ratio  in  the  horizontal  plane. 

An  investigation  of  the  effects  of  the  finite  length  of  the 
cylinder  shows  that  the  main  source  of  these  deep  nulls  is  the 
chosen  length  of  the  cylinder.  The  half-loop  radiator  and  various 
other  elements  were  analyzed  in  the  presence  of  an  infinite  cylinder 
by  the  methods  of  Chapter  III.  The  current  distribution  of  the 
radiating  element  is  first  calculated  by  the  moment  method  of 
Chapter  II.  This  segmented  distribution  is  then  used  to  weight  an 
infinitesimal  radiator  whose  fields  are  calculated  in  the  presence 
of  the  infinite  cylinder.  By  the  principle  of  superposition,  these 
fields  are  summed  to  give  the  far-field  pattern  of  the  radiator  on 
the  infinite  cylinder. 

Figure  4-8  shows  some  results  of  applying  this  technique  to 
the  half-loop  located  on  a  cylinder  of  infinite  length  with  diameter 
the  same  as  the  previous  model.  Comparing  these  patterns  with  Fig. 

4-7  shows  the  effects  of  the  finite  cylinder  length  on  the  theta 
component  of  far  electric  field  in  the  elevation  plane.  The  deep 
nulls  near  the  azimuthal  plane  have  vanished  over  the  frequency 
range,  but  the  desired  symmetry  about  the  horizon  is  only  slightly 
improved.  An  improvement  in  front-to-back  ratio,  especially  at  the 
lowest  frequency  is  also  apparent. 

To  relate  these  field  patterns  for  an  infinite  cylinder  to  the 
induced  voltage  patterns  of  the  input  ports  of  Fig.  4-3,  reciprocity 
may  be  employed.  The  theta  and  phi  components  of  the  far  electric 
field  for  a  half-loop  fed  from  bottom  to  top  are  calculated  and 
simply  multiplied  by  a  scale  factor  to  become  components  of  Port  1 
voltage.  Similarly,  fields  from  a  half-loop  fed  from  top  to  bottom 
become  components  of  Port  2  voltage.  These  port  voltages  may  be 
processed  as  described  previously  in  this  section  to  separate  the 
polarizations. 

Results  for  elevation  angles  of  0°,  10°,  20°  are  shown  in  Figs. 

4-9  through  4-11.  Comparison  of  these  results  near  an  infinite 
cylinder  with  those  of  Figs.  4-4  through  4-6  for  the  finite  length 
cylinder  demonstrate  the  effects  of  the  finite  cylinder.  The  theta 
components  of  the  induced  port  voltages  have  a  larger  front-to-back 
ratio  and  remain  more  constant  over  the  frequency  range  and  the 
range  of  elevation  angles  for  the  infinite  cylinder  case.  The  phi 
components  of  induced  port  voltages  are  only  slightly  more  stable 
for  the  infinite  cylinder  case. 

Thus,  this  section  has  demonstrated  the  method  of  separation  of 
polarization  of  the  incoming  signal  as  applied  to  a  half-loop  radiator. 
The  desired  DF  patterns  cf  Section  B  have  been  duplicated  over  approxi¬ 
mately  a  three-to-one  frequency  band  for  an  arbitrarily  polarized 
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Fig.  4-9(a).  Theta  component  of  received  voltage  pattern 

obtained  by  subtracting  voltage  at  Port  1  from 
Port  2.  Three  frequencies  are  shown.  The 
incident  wave  is  in  the  horizontal  plane.  The 
half-loop  elements  are  on  the  infinite  cylinder. 
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Fig.  4-9(b) .  Phi  component  of  received  voltage  pattern 
obtained  by  adding  voltage  at  Port  1  and 
Port  2.  Three  frequencies  are  shown.  The 
incident  wave  is  in  the  horizontal  plane. 
The  half-loop  elements  are  on  the  infinite 
cylinder. 
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Fig.  4-10(a).  Theta  component  of  received  voltage  pattern 

obtained  by  subtracting  voltage  at  Port  1  from 
Port  2.  Three  frequencies  are  shown.  The 
incident  wave  is  10°  above  horizontal  plane. 
The  half-loop  elements  are  on  the  infinite 
cylinder. 
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INDUCEO  PORT  VOLTAGES  VS-  PLANE  WAVE  INCIDENT  ANGLE 
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Fig.  4-1 0(b) .  Phi  component  of  received  voltage  pattern 
obtained  by  adding  voltage  at  Port  1  and 
Port  2.  Three  frequencies  are  shown.  The 
incident  wave  is  10°  above  horizontal  plane. 
The  half-loop  elements  are  on  the  infinite 
cylinder. 
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Fig.  4-1 1(b). 


Phi  component  of  received  voltage  pattern 
obtained  by  adding  voltage  at  Port  1  and 
Port  2.  Three  frequencies  are  shown.  The 
Incident  wave  is  20°  above  horizontal  plane. 
The  half- loop  elements  are  on  the  infinite 
cylinder. 


incident  signal  in  or  near  the  azimuthal  plane.  The  method  is 
presently  limited  to  small  elevation  angles  because  of  asymmetries 
in  the  elevation  plane  patterns  of  the  half-loop  radiator  in  the 
presence  of  the  cylinder.  Comparison  of  results  obtained  for  a 
finite  length  cylinder  with  those  of  an  infinite  length  cylinder 
shows  better  results  when  the  cylinder  is  electrically  long. 

F.  Input  Impedance  and 

fead i a t i on  Efficiency 

An  important  consideration  for  any  antenna,  transmitting  or 
receiving,  is  the  power  loss  involved  in  the  transfer  of  energy 
via  the  antenna.  Conventionally, this  loss  is  divided  into  two 
quantities.  One  is  mismatch  loss,  measured  as  VSWR,  due  to  the 
antenna  input  impedance  being  unequal  to  the  impedance  of  the 
connecting  line  or  cable.  The  second  is  radiation  efficiency  due 
to  losses  in  the  antenna  element  itself  such  as  omhic  losses.  For 
the  electrically  small  radiators  considered  in  this  study,  a  trade¬ 
off  between  antenna  efficiency  and  VSWR  was  clearly  established. 

A  direction  finding  system  would  be  most  concerned  with  antennas 
utilized  in  the  receiving  mode.  Radiation  efficiency  of  a  receiving 
antenna  can  usually  be  considerably  lower  than  for  a  transmitting 
antenna.  Extensive  investigation  of  antenna  input  impedance  versus 
antenna  efficiency  was  conducted  on  various  electrically  small 
radiators.  Perfectly-conducting  radiating  elements  tend  to  have 
highly  resonant  input  impedances.  That  is,  the  input  impedance  varies 
rapidly  over  a  frequency  band.  The  corresponding  VSWR  also  changes 
rapidly  and  can  be  quite  large,  even  over  a  small  frequency  range. 
Introducing  loss  into  the  radiator  in  the  form  of  reduced  conductivity 
or  lumped  loads  connected  with  perfectly-conducting  segments  damps  the 
rapidly  varying  input  impedance.  Thus,  radiators  formed  from  lossy 
wires  or  wires  which  are  loaded,  provide  a  means  for  stabilizing  the 
input  impedance  and  correspondingly  the  VSWR.  This  improvement  is 
made  at  the  expense  of  antenna  efficiency.  Resistive  losses  due  to 
the  decreased  conductivity  or  lumped  loads  greatly  reduce  the  antenna 
efficiency.  Various  loading  schemes  were  investigated  on  numerous 
radiating  elements;  however,  this  paper  will  present  results  for 
only  one  loading  scheme  as  applied  to  half-loop  radiators  of  Fig. 

4-3.  The  method  used  a  tapered  distribution  [7]  of  lumped  resistive 
loads  of  values  shown  in  the  figure.  These  values  are  not  picked 
to  optmize  any  particular  parameter,  but  to  demonstrate  the  possi¬ 
bilities  proper  loading  can  offer  in  terms  of  mismatch  versus  radia¬ 
tion  efficiency.  Figure  4-12  shows  a  Smith  Chart  mapping  the  input 
impedance  over  a  four-to-one  frequency  band  when  normalized  to  300  n. 
The  VSWR  remains  below  1.65.  Table  4-1  shows  the  corresponding 
antenna  efficiencies  for  this  frequency  range. 
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CHAPTER  V 


SUMMARY  AMO  CONCLUSIONS 

To  summarize,  a  method  of  direction  finding  has  been  proposed 
and  demonstrated.  The  method  evolved  in  the  process  of  surmounting 
the  constraints  of  the  problem.  The  antenna  radiating  elements 
were  to  be  electrically  small  and  located  on  a  cylinder  whose 
electrical  dimensions  were  larger.  This  caused  the  component  field 
pattern  shapes  and  directivity  to  be  limited  to  two  basic  kinds: 
the  figure  eight  phi  pattern  and  the  broad  single  lobed  theta 
pattern.  When  considered  independently,  these  patterns  were  shown 
to  be  capable  of  direction  finding;  however,  practical  elements 
investigated  induced  both  simultaneously.  By  correctly  combining 
properly  oriented  radiators,  the  independent  patterns  were  re¬ 
covered  for  arbitrarily  polarized  DF  signals  in  or  near  the 
azimuthal  plane. 

Further,  the  system  was  to  perform  over  a  broad  frequency  range 
with  the  frequency  of  the  incident  signal  unknown.  Processing  of  the 
incoming  signals  was  to  be  of  the  most  elementary  kind,  that  is,  a 
simple  processing  algorithm  was  to  work  for  all  signals  in  the 
operating  band.  These  constraints  required  that  the  pattern  shapes, 
including  the  location  of  maxima  and  minima,  remained  constant  over 
the  frequency  band.  These  necessary  requirements  were  achieved  and 
demonstrated  for  a  particular  radiating  element  over  a  three-to-one 
frequency  band. 

An  antenna  system  operating  over  a  frequency  range  has  the  added 
difficulty  of  maintaining  a  reasonably  well-behaved  input  impedance 
so  that  mismatch  losses  are  minimized.  Perfectly-conducting  radiators 
of  the  type  investigated  have  rapidly  varying  input  impedances  over 
the  frequency  band.  A  method  of  tapered  loading  has  been  demonstrated 
which  effectively  damps  this  variation  at  the  sacrifice  of  antenna 
efficiency.  Thus,  a  method  using  electrically  small  radiators 
properly  formed  into  an  array  around  a  cylinder  has  been  demon¬ 
strated  to  be  capable  of  direction  finding  over  a  frequency  band 
with  minimal  signal  processing,  even  when  the  frequency,  phase  and 
polarization  of  the  incoming  signal  are  unknown. 

The  investigation  of  the  DF  system  was  facilitated  by  the  use 
of  modern  numerical  techniques  in  antenna  analysis.  Radiating 
elements  in  the  presence  of  the  cylinder  were  analyzed  using  two 
approximate  techniques.  One  was  a  specialized  moment  method  in 


which  the  antenna  structure  was  wire-grid  modeled,  end  piecewise- 
sinusoidal  expansion  and  weight  functions  were  assumed  and  reaction 
enforced.  This  method  was  applied  to  radiators  on  finite  cylinders. 
The  second  method  was  a  hybrid  technique  which  used  currents  calcu¬ 
lated  by  the  first  method  to  weight  infinitesimal  elements  near  an 
infinite  cylinder  whose  far  fields  were  found  by  modal  solutions. 
Integration  over  the  element  apertures  gave  the  total  fields.  The 
accuracy  of  the  wire-grid  modeling  method  has  been  well  documented 
in  the  literature.  The  second  method  uses  exact  equations  for  the 
fields,  the  approximation  being  the  weighting  factor  calculated  from 
a  discrete  current  distribution  generated  on  a  finite  cylinder.  The 
accuracy  of  this  method  is  very  good  when  the  elements  are  loaded 
and  sufficient  current  samples  are  used  over  the  element. 

Several  limitations  of  the  specific  radiator  chosen  to  demon¬ 
strate  the  DF  system  concept  were  discussed.  The  half-loop's  phi 
component  of  electric  field  in  horizontal  plane  had  a  maximum  which 
drifted  slightly  over  the  frequency  range.  The  theta  component  of 
electric  field  in  the  elevation  plane  was  asymmetric  about  the 
horizon,  causing  the  polarization  separation  technique  to  work  for 
only  small  angles  of  elevation.  These  limitations  for  this  particular 
radiator  do  not  preclude  the  use  of  this  concept  for  systems  which 
need  a  broader  bandwidth  or  higher  elevation  angles.  A  radiating 
element  would  need  to  be  designed  with  more  symmetric  elevation  plane 
patterns  and  less  pattern  shift  over  a  larger  frequency  bandwidth. 

Results  of  the  investigations  of  small  radiators  on  the  conduc¬ 
ting  cylinder  suggest  several  topics  for  further  research.  First, 
a  search  for  a  practical  radiating  element  which  is  sensitive  to  phi 
polarized  signals' and  blind  to  theta  components  should  be  completed. 
One  possible  candidate  is  a  small  loop  located  in  the  horizontal 
plane.  Elements  with  more  symmetric  elevation  plane  patterns  could 
be  utilized  to  extend  the  method  for  higher  elevation  angles.  Design 
curves  could  be  calculated  so  that  for  given  loss  specifications  an 
optimum  mismatch  versus  efficiency  trade-off  could  be  obtained  for 
various  loading  schemes.  Finally,  experimental  verification  of  the 
computer  calculated  results  should  be  carried  out.  This  would  be 
especially  desirable  and  satisfying  for  a  system  whose  specifications 
are  fixed,  including  frequency  bandwidth,  tolerable  power  loss, 
necessary  accuracy,  and  conformal  geometry  constraints. 
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